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ABSTRACT : Two important ingredients of nuclear shell-structure, kinetic energy and 

> : 

spin-orbit splitting, are studied as a function of orbital angular momenta £ and binding 
energies, when binding energies of neutrons decrease towards zero. If we use the standard 
parameters of the Woods-Saxon potential in (3 stable nuclei and approach the limit of 
zero binding energy from 10 MeV, the spin-orbit splitting for n—1 orbitals decreases 
considerably for £=1, while for £>2 little decreasing is observed in the limit. In contrast, 
the kinetic energy decreases considerably for £ < 3. The smaller the £ values of orbitals, the 
larger the decreasing rate of both kinetic energy and spin-orbit splitting. The dependence 
of the above observation on the diffuseness of potentials is studied. 



PACS numbers : 21.10.Pc 



1 Introduction 



One-particle shell structure is well studied for nuclei along the (3 stability line. Using recent 
radioactive nuclear ion beam facilities it has become possible to study the structure of 
exotic nuclei away from the /3 stability line, in which least-bound nucleons have much 
smaller separation energies compared with those in f3 stable nuclei. Some interesting 
change of shell structure is expected and/or observed, when we move away from f3 stable 
nuclei to drip line nuclei. This change is made due to the presence of loosely-bound 
nucleons as well as exotic values of isospin for a given mass number. For example, available 
experimental information on the nucleus l^^-Szoi which is a singly-closed-shell nucleus 
in the traditional terminology, seems to indicate that the nucleus is deformed |], @. 
Furthermore, the creation of a new magic number N=16 near the neutron drip line (T 2 >3) 
is recently suggested || . In the present work we study the two basic ingredients of nuclear 
shell-structure, kinetic energy and spin-orbit splitting, for relatively simple potentials, 
having the spirit similar to that of ref. Q. In available publications certain aspects of 
shell structure unique in drip line nuclei have been studied using self-consistent Hartree- 
Fock (HF) calculations of various types. 

Using Skyrme-type HF calculations, in ref. we have studied neutron one-particle 
level scheme with a fixed mass number as a function of neutron number, including one- 
particle resonant levels. It is seen that calculated one-particle levels with lower i are 
systematically pushed down relative to those with higher £, as we approach the neutron 
drip line. On the other hand, in ref. || the reduction of the spin-orbit potential in light 
neutron rich nuclei, especially the isospin-dependence of the spin-orbit term is discussed, 
using the relativistic mean field theory. In the Skyrme-type HF calculations of ref. |7j 
it is found that the estimated IP1/2-I-P3/2 spin-orbit splitting of the light mirror nuclei, 
gC 2 and 2He 6 , is about 30 percent smaller than the one in the (3 stable nucleus g 2 C 6 . The 
difference of the isospin-dependence of spin-orbit term in the relativistic and traditional 
Skyrme-type HF calculations is certainly important || to be examined. 



A systematic study of both kinetic energy and spin-orbit splitting as a function of 
one-particle binding energy, orbital angular momentum and the shape of the potential 
is essential in understanding the variation of the shell structure as we approach the drip 
lines. In self-consistent HF calculations not only the orbitals with less binding energies 
are occupied but also the effective diffuseness of one-body potentials is changed due to 
the occupation, as we approach the drip lines. In the following the variation of shell 
structure is studied for potentials with a given diffuseness. We note that in usual shell 
model calculations harmonic oscillator wave functions are used and single-particle energies 
are input parameters. In order to give appropriate input parameters for nuclei away from 
the (3 stability line, it is important to know characteristic features of both the variation 
of kinetic energy and that of the strength of spin-orbit potential when approaching the 
drip lines. In ref. || the isospin-dependence of kinetic energies in light neutron-rich 
nuclei is suggested to be used in shell-model calculations. We will later comment on this 
suggestion. In the present work we limit ourselves to the study of lighter nuclei, since at 
the moment experimental data on nuclei close to the neutron drip line may be obtained 
only for lighter nuclei. 

In Fig. [l]a we plot one-particle energy eigenvalues for neutrons in the Woods-Saxon 
potential 



choosing the lp, 2s, Id, 2p and If orbitals. In the expression (Q) A c expresses the nucleon 
Compton wave-length. Used parameters are 

R = r A 1/3 with r = 1.27 fm 




a 



(1) 



together with the spin-orbit potential 





a 



0.67 fm 



(3) 



and N=Z in the depth of the Woods-Saxon potential 



A 

and in the strength of the spin-orbit potential |10 



N - 

U = (50 - 32 — ) MeV (4) 



N-Z\ 

v = 32 - 20 — (5) 
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The energies of neutron orbitals calculated using almost the same parameters as in Fig. 
^a are found for A>20 in p. 239 of ref. |10j. Due to the presence of less-bound neutrons the 
diffuseness of the potential may become larger as we approach neutron drip lines. However, 
just using a larger diffuseness may not always correspond to the realistic situation. This is 
partly because the major part of the total potential comes from the core particles, which 
are well bound and are the majority. Furthermore, it is partly because the proton core, 
which makes a more important contribution to the neutron potential than the neutron 
core, is more deeply bound in nuclei closer to neutron drip lines. As a numerical example 
with a larger diffuseness, in Fig. |I]b we show the one-particle energy eigenvalues for the 
Woods-Saxon potential with the same parameters as those used in Fig. |l]a except for 
doubling the diffuseness, a=1.34 fm. 

In Figs. |l]a and |T]b we observe : (a) the A-dependence of the eigenvalues E ni in the 
energy region of E n g> — 15 MeV depends very much on the orbital angular momentum ; 
(b) as is seen, for example, from the crossing of the eigenvalues of lds/2 and 2si/2 orbitals 
at £ n ^=-0.71 MeV in Fig. 0a and £^=-7.21 MeV in Fig. 03, it is known that the 
slope of eigenvalues for lower £ orbitals becomes relatively much less, as the eigenvalues 
approach zero ; (c) The spin-orbit splitting with a given (n, £) may generally be expected 
to become larger (~A -1 / 3 ) as the mass number decreases. In contrast, the spin-orbit 
splitting between lpi/2 and IP3/2 (or 2pi/ 2 and 2p 3 / 2 ) does not increase but it rather 
decreases when E ni goes up from —5 to MeV. In Fig. |l|b the decrease is observed also 
for the spin-orbit splitting of Id orbitals; (d) For the potential with a larger diffuseness 
the features mentioned above, (a)-(c), are more pronounced. We note that the prominent 



difference of one-particle shell structure, such as the level order and the energy gap, 
between Figs. 0a and |l]b is observed especially for E ne > —10 MeV. 

In ref. [ll]| we studied the coupling of particles near threshold to shape oscillations 



and showed the evaluated matrix elements of 

VAr) = *^ (6) 

where U(r) expresses the one-particle potential while R represents the nuclear radius. 
Noting that the kinetic energy of the spherical system in the absence of spin-orbit potential 
is equal to the expectation value of 

rdU(r)_ 

2 dr U 
due to the virial theorem, for the square-well potential the expectation value of (0) is 
equal to twice the kinetic energy. Furthermore, the expectation value of 

i<m (8) 

r dr 

which is the radial dependence of the spin-orbit potential, is equal to 1/R 2 times that 
of <M) in the case of the square-well potential. Thus, for the square-well potential Fig. 
3a of ref. [ I l]] already contains all information on both the kinetic energy and the spin- 
orbit splitting. In the present work we investigate those two quantities for Woods-Saxon 
potentials compared with the square-well potential, studying the dependence on energy 
especially for E n ^=— 10^0 MeV and on angular momenta. Though it can be disputed 
whether the radial derivative of the density or that of the potential should be used for 
the radial dependence of the spin-orbit potential, for simplicity we use the form in 

In sect. 2 we present numerical results and discussions, while conclusions are given in 
sect. 3. 



2 Numerical results and discussions 



In the present section we do not include the spin-orbit potential in the one-body potential 
for simplicity. We evaluate the expectation value of the operators, 

r dU(r) 



Uq dr 

and 



(9) 



1 iV(r) Q) 



r dr 

for neutron orbitals in Woods-Saxon potentials. The expectation values of (^) is equal 



to those plotted in Fig. |3|a of JTTJ in the case of the square-well potential. Note that the 
quantity on the x-axis in Fig. [I]a of ref. |TT[ for the square-well potential has the same 
meaning as that on the x-axis in Fig. and of the present paper, though the notation 



used for the square- well potential in ref. ]TTJ] was different. The expectation value of (JD 
is proportional to the kinetic energy, while that of (|10D to the spin-orbit splitting for a 
given orbital angular momentum i treating the spin-orbit potential by perturbation. We 
examine those expectation values for a given (n,£) value as a function of E n i/Uo where 
E n £ expresses one-particle energy eigenvalues. 

In Figs. 0a and 0b the expectation value of the operator (0) is plotted for the Woods- 
Saxon potential with the diffuseness a=0.67 and 1.34 fm, respectively, while in Fig. 0c 
we show the same quantity for the square-well potential, 

{—U Q for r < R 

(11) 
for r > R 

a part of which was taken from Fig. 0a of ref. [11]. In the limit of infinite square- well 
potential (Uo^oo), the kinetic energy is equal to Uo+E n £ and, thus, the expectation value 
of (0) multiplied by U is given by 2(XJ +E n i) for all orbitals with various (n,£) values. 
The expectation value for the infinite square-well potential is shown in Fig. 0c by the thin 
straight line. Since the relation between the quantities in the x- and y-axis multiplied 
by U is expressed by a straight line, dividing both quantities by any given "U " does 



not change the straight line. Thus, we plot the line for the infinite square-well potential 
also in Fig. |2|c. Note that one-particle wave-functions cannot extend to outside of the 
potential with Uo— >oo and, thus, the reduction in the increasing rate of kinetic energy 
does not occur as (l+(E n ^/Uo)) increases. 

It should be noted that along the curve with a given (n,£) the mass number mono- 
tonically decreases as (l+(E n i/XJ )) increases. For a very large value of U , one may 
approximately write (see eq. (|T3|)) 

Uo + E ni= C f 2 (12) 

where c n ^>0, U >0 and E n £<0, while R expresses the radius of the system. Thus, in 
order to obtain larger eigenvalues E n g, namely larger values of (l+(E n £/Uo)), keeping a 
fixed value of Uo , one must decrease values of R (or the mass number A). 

The vertical dotted line of Fig. 2 at l+(E n€ /U )=0.8, namely E ne =-10 MeV for U =50 
MeV, indicates the approximate Fermi level in (5 stable nuclei. When the eigenvalue 
E n £— >-0 exists for the potential with a finite range, the expectation value of @ for (n, £=0) 
orbitals in the finite-well potential approaches zero, since due to the absence of centrifugal 
barrier the £=0 neutron wave functions can extend up to infinity. For the square-well 
potential the expectation value of @ at E n ^=0 is independent of n and is equal to 2/3, 
6/5 and 10/7, for the £=l, 2 and 3 orbitals, respectively [fl"T|. In contrast, the expectation 



value for Woods-Saxon potentials depends on radial node n and is almost a factor two 
(four) smaller for the potential in Fig. (0b). It is also noted that for more diffuse 
potentials the absolute magnitude of kinetic energy of a given (n, £) orbital is considerably 
smaller for almost all values of E n g. 

The effect of the diffuse surface of the potential on kinetic energy is seen from the 
comparison of the curves in Figs. |2|a and |2|b with those in Fig. 0c for the finite and 
infinite square-well potential. The considerable deviation from the straight line occurs 
first at smaller values of | E n g | for less diffuse potential and for larger £ orbitals. In other 
words, the kinetic energy starts to decrease already at relatively larger binding energy 



for more diffuse potentials. This is because for more diffuse potentials the wave functions 
can more easily extend to the outside of the potentials. The extension is easier for the 
orbitals with smaller orbital angular momenta because of the lower centrifugal barriers. 
As seen in Fig. |^c, for the square-well potential the kinetic energy starts to decrease first 
at very small binding energies, even for £=0 orbitals. 

The deviation of kinetic energy of (n, £) orbitals from the monotonic increase as E n ^ 
increases from —10 to MeV depends somewhat on the radial node n. The larger the 
quantum number n, the larger the variation of the kinetic energy for E n ^=— 10 — > MeV. 
This is understood because at (1+E n ^/U )=0.8, namely E n £=— 10 MeV, the kinetic energy 
of orbitals with larger n is larger than that with smaller n, while at E n ^=0 the kinetic 
energy of orbitals with a given I is almost independent of n. The latter is independent of 
n for the square-well potential. 

In ref. the isospin dependence of kinetic energies in very light neutron-rich nuclei 
is suggested, in which the isospin dependence for the 2si/ 2 and ld 3 / 2 state is taken to be 
the same as that for the ld.5/2 state. The origin of the isospin-dependence is certainly 
traced back to the decreasing binding energy of least bound neutrons as well as the 
possible larger surface diffuseness of the potential, when the isospin increases in those 
light neutron-rich nuclei. However, from our study shown in Fig. ^| we note that the 
binding-energy dependence of kinetic energy for the 2s and Id states are very different 
for a given potential, especially for E n ^=— 10 — > MeV. Observe the difference of the 
behaviour between the solid and dashed curve in Fig. Q for (1+E n ^/Uo)=0.8— >1.0 . 

In Figs. ||a and ||b the expectation value of the operator (|i0|) is plotted for the Woods- 
Saxon potential with the diffuseness a=0.67 and 1.34 fm, respectively, while in Fig. |3]c 
we show the same quantity for the square-well potential. For the square-well potential 
the expectation value of ([H]) depends on some extra parameter to specify, either Uo or 
R, other than the only parameters unique in the model, Uo-R 2 and E n iR 2 [0. We have 
used U =50 MeV. 



The thin curves in Fig. ||]c are estimated for the infinite square-well potential. For the 
(U — >oo) square-well potential the kinetic energy is equal to 

h 2 x 2 

U + E ne = (13) 

where x n g expresses the n-th zero point of the spherical Bessel function je(z). Thus, the 
expectation value of the operator in ( |10"D multiplied by Uo is written as 

2{U + E nl ) 4M(C/ + E nl f 



R 2 h 2 x 2 



(14) 



Variables in the x- and y-axis in Fig. 3c are those in ( |T3| ) and (|TJ|) divided by "Uo", 
however, the division by "Uo"=oo does not make sense. Since the relation between the 
quantities in (|13|) and (|14"D is not linear, the shape of the curve depends on the "U " value 
used for evaluating variables in the x- and y-axis. Since our purpose is to compare the 
result of finite potentials with that of the infinite potential, we choose this "Uo" value so 
that at (l+(E n ^/Uo))=0.8 the thin curve coincides with the corresponding curve for the 
finite square-well potential. It is seen that only in the region (l+(E n f/U )) > 0.8 the curve 
for the finite square- well potential differs appreciably from that for the infinite square-well 
potential. 

The plotted quantity in Fig. |3] times (2£+l) is proportional to the spin-orbit splitting. 
For a given E n e value the spin-orbit splitting for the lp orbital is much larger than that for 
the 2p orbital, since the former belongs to lighter nuclei. For the purpose of comparison 
of expectation values of fllOD , we plotted calculated values also for 2s orbitals, though no 
spin-orbit potential acts on s-orbitals. 

The diffuse surface effect of potentials on spin-orbit splitting is seen from the compari- 
son of curves in Figs. 3a and 3b with those in Fig. 3c for the finite and infinite square-well 
potential. Due to the different r-dependence, the effect of the diffuse surface on kinetic 
energy starts to appear already at deeper bound orbitals compared with spin-orbit split- 
ting. As already exhibited in Fig. |I|b, for the Woods-Saxon potential with a=1.34 fm 
the spin-orbit splitting of Id orbitals decreases also appreciably as E n £ changes from —10 



to MeV. The decrease corresponds to the variation of the dashed curve in Fig. §b for 
(l+E n ,U )=0.8^1.0 . 

The variation of spin-orbit splitting of orbitals with a given £ for E n £=— 10— >0 MeV 
depends moderately on the radial node n. We also note that the comparison of spin-orbit 
splittings of the orbitals with a given £ but with different n means the comparison made 
for nuclei with different mass number. 

In order to simulate the spin-orbit potential in which the radial derivative of the den- 
sity is used instead of that of the potential, we have also made the numerical calculation 
using r =(0.9)(1.27)=1.14 fm in U(r) of (|TUD while keeping one-particle wave-functions un- 



changed. The dependence of the expectation value of ([LT]) on either £ or n for E ni =— 10^0 
MeV is found to be almost the same as what is described above. 

3 Conclusions 

In order to study the variation of one-particle shell-structure due to the surface diffuseness 
of nuclear potential, the radial matrix-elements of both kinetic energy and spin-orbit 
potential of Woods-Saxon potentials are estimated and compared with those of the infinite 
and finite square-well potential. We are especially interested in the energy region where 
binding energies of neutrons decrease from 10 to MeV, having nuclei close to neutron drip 
lines in mind. In the infinite square-well potential neither kinetic energy nor spin-orbit 
splitting ever decreases as a function of E n £. For the finite square-well potential with 
Uo=50 MeV a decrease is observed for E n £=— 10^0 MeV, depending on orbitals. For 
orbitals with smaller £ values both quantities decrease more strongly as binding energies 
approach zero. The decrease is stronger for more diffuse potentials. Due to the different 
r-dependence, the effect of the diffuse surface on kinetic energy starts to appear already 
at deeper bound orbitals compared with spin-orbit splitting. 

In traditional shell model calculations carried out for neutron drip line nuclei, it is 
important to take into account the variation of both kinetic energy and spin-orbit split- 



ting in the input one-particle energies, which is studied in the present work. Not only the 
binding energy dependence of those quantities but also the I dependence is very important. 
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Figure 1: Energies of neutron orbits : (a) for the Woods-Saxon potential with parameters, 
©j (H)> (H) an d N=Z ; (b) the same as (a) except for a=1.34 fm. 
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Figure 2: The expectation values of the operator (|5|), which is proportional to kinetic en- 
ergy of respective one-particle orbitals, as a function of (l+(E n ^/Uo)) where E ni expresses 
the one-particle energy eigenvalue of the orbital (n, €). For U =50 MeV (l+(E n ^/U ))=0.8 
and 1.0 mean E n e=— 10 and MeV, respectively, (a) for the Woods-Saxon potential with 
parameters, (§), @, v=0 and N=Z ; (b) the same as (a) except for a— 1.34 fm ; (c) for 
the square- well potential (PH). By the thin straight line in Fig. 0c the expectation value 
of the operator (§) is exhibited for the infinite square- well potential. The thin line is the 
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Figure 3: The expectation values of the operator fliTf), which in the perturbation treatment 
of spin-orbit potential is proportional to the spin-orbit splitting of the £ orbital when 
multiplied by (2m), as a function of (1+(E^/U )). For U =50 MeV (l+(E n€ /U ))=0.8 
and 1.0 mean E n e=— 10 and MeV, respectively, (a) for the Woods-Saxon potential with 
parameters, @), v=0 and N=Z ; (b) the same as (a) except for a=1.34 fm ; (c) for the 
square-well potential ([0]) with Uo=50 MeV. By the thin curves in Fig. |^c the expectation 
value in the infinite square-well potential is shown for the lp and Id orbitals. See the text 

for Hptnilc: 



